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We investigate the late-time tails of self-interacting (massive) scalar fields in the spacetime of dila- 
ton black hole. Following the no hair theorem we examine the mechanism by which self-interacting 
scalar hair decay. We revealed that the intermediate asymptotic behavior of the considered field 
perturbations is dominated by an oscillatory inverse power-law decaying tail. The numerical sim- 
ulations showed that at the very late-time massive self-interacting scalar hair decayed slower than 
any power law. 
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The late-time evolution of various fields outside a collapsing body plays an important role in major aspects of black 

hole physics, as the no-hair theorem y| and the mass inflation scenario Wheeler coined the metaphoric dictum 

O" 1 black holes have no hair which means that regardless of the specific details of the collapse or structure and properties 

of the collapsing body a stationary black hole characterized by its mass, charge and angular momentum emerged in 

, the resultant process. Therefore it is interesting to study the mechanism responsible for the decay of black hole hair. 
• i— I . 

The neutral external perturbations were first studied in || and it was found that the late-time behavior is dominated 
by the factor £~( 2 '+ 3 ) 5 for each multipole moment I. The neutral perturbations along null infinity and along the future 
event horizon were studied in Ref. Their decays were governed by the power laws u _ ^ +2 - ) and v~( l+3 \ where u 
and v were the outgoing Eddington-Filkenstein (ED) and ingoing ED coordinates. Bicak Q considered the scalar 
perturbations on Reissner-Nordtrom (RN) background and found that for | Q |< M the relation £~( 2 '+ 2 ) ; while for 
| Q |= M the late-time behavior at fixed r was governed by t~( l+2 \ 

Hod and Piran ||- |8) studied the late time behavior of a charged massless scalar field in RN spacetime. Among 
all the conclusion was that a charged hair decayed slower than a neutral one. 

The problem of late-time tails in gravitational collapse of a self-interacting (SI) fields in the background of 
Schwarzschild solution was reported by Burko |)j and in RN spacetime, at intermediate late-time, was considered 
in Ref. pj. At intermediate late-time for small mass m the decay was dominated by the oscillatory inverse power 
tails £ - ('+ 3 / 2 ) sin(mi). This analytic prediction was verified at intermediate times, where mM < mt < l/(mM) 2 . In 
Ref. |pd|] the nearly extreme RN spacetime was considered and it was found analytically that the inverse power-law 
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behavior of the dominant asymptotic tail is of the form i~ 5 / 6 sin(mi), independent of I. The asymptotic tail behavior 
of SI scalar field was also studied in Schwarzschild spacetime [Q. The oscillatory tail of scalar field has the decay 
rate of i -5 / 6 at asymptotically late time. 

Nowadays, it seems that the superstring theories are the most promising candidates for a consistent quantum theory 
of gravitation. In Ref. |L3| the asymptotic evolution of a massless charged scalar field in the background of dilaton 
black hole was studied. Using both numerical and analytical methods the inverse power-law relaxation was revealed. 
It turned out that the charged hair decayed slower than a neutral one. 

In our work we shall generalize our previous considerations [ p"3| and discuss the SI scalar field behavior in the 
spacetime of dilaton black hole. In Sec. II we gave the analytic arguments concerning the intermediate and the very 
late-time behavior of SI scalar field in the background of the considered black hole. Then, in Sec. Ill we treated the 
problem numerically. Sec. IV will be devoted to a summary and conclusions. 

II. THE EINSTEIN-MAXWELL-DILATON EQUATIONS 

In this section we shall study the evolution of SI (massive) scalar field ip around a fixed background of electrically 
charged dilaton black hole. The wave equation for the field is given by 

V^VV - m 2 ^ 2 = 0. (1) 

where m is assumed to be real. 

The metric of the external gravitational field will be given by the static, spherically symmetric solution of Eqs. of 
motion derived from the low-energy string action (see e.g. [ljj]). The action has the form as follows: 

S = J d*x [R - 2(V0) 2 - e~ 2 *F 2 ] , (2) 

where <f> is the dilatonic field and F a p — 2V[ a Ag]. The metric of electrically charged dilaton black hole G3] implies 



ds 2 = -(l- —) dt 2 + -. — — +r(r-%\ (d9 2 + sin<?# 2 ). (3) 
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The event horizon is located at r+ = 2M. For the case of r_ = ^ we have another singularity but it can be ignored 
because of the fact that r_ < r + . The dilaton field is given by e 2 ^ — e~~ 2 ^ (l — — ), where <j>o is the dilaton's value 
at r — > oo. The mass M and the charge Q are related by the relation Q 2 = r+ ^~ e 2 ^ . 
Defining the tortoise coordinate y, as dy = ^ dr ^ M ^ one can rewrite (||) in the form 



For the spherical background each of the multipole of perturbation field evolves separetly so for the scalar field in the 
form ip — J^i m V'LtA r )^j m (^; 4>)/R( r )i one nas the following equations of motion for each multipole moment 

ttt ~ ^.,w + Vip = 0, (5) 

where 
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By i? we denoted R = \Jr yr — ^-J and ' is the derivative with respect to the r-coordinate. 

In order to analyze the time evolution of SI held in the background of dilaton black hole we shall use the spectral 
decomposition method jL5| . The time evolution of SI scalar field is given by 



^(y,t) = / dy' 



G(y, y'\ t)My', 0) + G t (y, y'; t)<P(y' , 0) 



for t > 0, where the Green's function G(y, y'; t) is defined by 
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G(y,y';t)=S(t)S(y-y'). 



(7) 



(8) 



dt 2 dy 2 

Our main task will be to find the black hole Green function. The first step in finding G(y, y'; t) consists of reducing 
(^J) to an ordinary differential equation. To do it one can use the Fourier transform 

(9) 



(10) 



G(y,y';cj) = / dtG(y,y';t)e^. 
Jo- 

This transform is well defined if Im lo > 0, while the corresponding inverse transform has the form 



G(y,y';t) = 
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duG{y,y';uj)e- 
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for some positive number e. 

The Fourier's component of the Green's function G(y, y'; to) can be expressed in terms of two linearly independent 
solutions for homogeneous equation, namely 

(j^ + ^ - V)^i = 0, i=l,2. (11) 

The boundary conditions for %pi are described by purely ingoing waves crossing the outer horizon H+ of the dilaton 
black hole ipi ~ e - lu] v as y — > — oo while ip2 should be damped expotentially at ?'+, namely ip2 — eT ~ m y at 
y -> oo. 

In order to find ipi we consider the wave Eq.([ll]) of SI scalar field and transform it in such a way that the 
Coulomb and Newtonian (1/r) potentials will dominate at large distances. One can introduce an auxiliary variable 
£ in such a way that £ = A(r)tj), where A(r) 2 = 1 — 2M. So in terms of the variable £ Eq.(|ll]) can be rewritten as 
follows: 

'1(1 + 1) R" . 2MR' 



d £ A rr . 1 

S ,rr_ £ _^ 

rfr 2 A A* 



u 2 - A 2 1 + ^A(r) 2 
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e = o. 



(12) 



Now we expand Eq.(|l2|) in power series of and ®, neglecting terms of order O 
following expression: 



dr 2 



, , 4Mw 2 - 2Mto 2 1(1 + 1) 
lj 2 - m 2 H v n ' 



£ = 0. 



Thus we arrive at the 



(13) 



If one further assumes that the observer and the initial data are in the region where r < ^J^yi and one shall be 



interested in the intermediate asymptotic behavior of SI scalar field r < t < 



M 



(Mr, 



then we get 
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dr 2 



2 2 
ijj — m 



1(1 + 1) 



e = o. (M) 



The above approximation neglects the backscattering of the SI field from asymptotically far regions. As in Ref. || in 
our case of the scalar field perturbations on dilaton black hole wc get the dependence on the field's parameter (the 
mass of the field). These perturbations do not depend on the spacetime parameters (Q and M). 

The procedure of getting the solution of Eq. ( |l4| ) is the same as described in || so we refer the reader to this work. 
We only write the final form of the ip dependence. The intermediate asymptotic behavior of the SI field at a fixed 
radius has the form ^ 

ijj ~ S i n mt, (15) 

while the intermediate behavior of the considered field at dilaton black hole outer horizon H + is dominated by 

ip ~ v ~( l +i) sxnmt. (16) 

They are dominated by an oscillatory power law tails. 

Recently the late-time behavior of the SI field in Schwarzschild background was studied in Ref. ]l^| . Following 
their footsteps we shall write the wave Eq. ( |l4| ) using the non-dimensional variable x = r/2M . Then we reach to the 
analytically tractable equation of the form 



d 2 tp 1 dtp 

dx 2 x(x — 1) dx 



AM 2 uj 2 x 2 1(1 + 1) 1 Am 2 M 2 x 



r + (>[^) 0. (17) 



x — l) 2 x(x — 1) x 2 (x — 1) x — 1 

The procedure of getting the solution to Eq.(|l7|) was presented in details in Ref. Jj^ |. Here we quote the main 
conclusion, namely the smaller value of raM is the later the i -5 / 6 tail begins to dominate. It is true for the range 
of parameter rat ^> jj^jjyi- For rat 3> raM it was found that the larger value of raM is the later the i~ 5 / 6 tail 
begins to dominate. Of course these conclusion are valid for small Q. In the next section we investigate this problem 
numerically. 

III. NUMERICAL RESULTS 

We numerically analyzed Eq. (^) using method described in We transformed Eq. (||) into (u, v) coordinates 

4^ l1M , + Vi/> = 0, (18) 

and solved it on uniformly spaced grid using explicit difference scheme. As was pointed out previously the late time 
evolution of a massive field is independent of the form of the initial data. To start our calculation we used a Gaussian 
pulse of the form as follows: 

i,(u = Q,v) = Aexp(- ( V ~^ ) (19) 



Linearity of the Eq.(||) leaves us a freedom to choose the amplitude A, and for our purposes we used A = 1. The 
rest of the initial field profile parameters were vq = 50 and a — 2. The black hole mass and charge are set equal to 
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M = 0.5 and Q = 0.45, respectively, and mass of the field is m = 0.01. We have studied the behavior of the field ip 
on two hypersurfaces: the future timelike infinity i + (approximated in our calculations by the field at fixed radius 
y = 50), and the black hole future horizon H + (approximated by the field on the null surface u = 10 4 ). First, we have 
examined I = modes. The results are presented in Fig. [l| After initial period dominated by prompt contribution and 
quasi- normal ringing the field establishes a definite oscillatory pattern. The amplitudes of the field decrease according 
to power-law with index —1.51 for both i + (top curve) and H + . These exponents were calculated by fitting a straight 
line to the oscillations extrema within the time range 2 x 10 3 -^7 x 10 3 to avoid influence of the initial perturbation and 
the late time tails. Analytically estimated from Eqs . (|l5|)- (|l6|) value for the index is —1.5 providing excellent agreement 
with numerical results. The period of the oscillations (the distance between two neighbor maxima) is 314 ±0.5 which 
means less than 0.2% discrepancy from analytically expected period T = n/m ~ 314.5 
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FIG. 1. Evolution of the field | tjj | for I = 0, M = 0.5, Q = 0.45 and m = 0.01. Upper curve represents the field at fixed 
radius y — 50 (approximation of the future timelike infinity i+) as a function of t. Lower curve represents the field at u — 10 4 
(approximation of the black hole future horizon H+) as a function of v. The power-law exponents are —1.51 for both i+ and 
H+. The period of the oscillations is T = n/m ~ 314.5 to within 0.2%. 

Secondly, we considered the dependence of the SI scalar field on the multipole index I. The results are presented in 
Fig. ||. The behavior of the field on the future timelike infinity is calculated for I — 0,1,2 and 3. The numerical values of 
the power-law exponents are —1.51, —2.52, —3.53, and —4.51, respectively. According to Eq.([l5|) these exponents equal 
to —1.5, —2.5, —3.5, and —4.5, thus again the agreement between numerical calculations and analytical predictions 
is excellent. The oscillations period is T = n/m ~ 315 to within 0.5%, in agreement with the predicted value. 
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FIG. 2. Evolution of the field \ip\ on the future timelike infinity i+ (approximated by ip(y = 50, t)) for different multipoles 
I = 0, 1, 2 and 3. The power-law exponents are —1.51, —2.52, —3.53, and —4.51, respectively. The period of the oscillations is 
314.5 to within 0.5% (for the worst case of I = 3). 

Finally, we have checked the late-time behavior of the field. In Fig. || we present the field on the future timelike 
infinity i + and on the black hole horizon H + as a function of time for I = mode (for presentation purposes we plot 
only maxima of the oscillations) . As in Jl0| our numerical scheme reveals the fact that after the intermediate phase of 
the power-law decay the amplitude of the field decreases slower than any power-law. In the very late approximation 
the change of the tail will be result of the dominant backscattering due to the spacetime curvature. This effect is 
beyond the flat space approximation as was remarked in Jll[ . 
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FIG. 3. Late-time behavior of the field | ip \ m ax on dilaton black hole. The mass of the field m = 0.05, the upper curve 
represents the field at future timelike infinity (y — 50) as a function of time t. The lower curve along the black hole horizon as 
a function of v. The period of oscillations is 63.0 to within 0.8%. The thin solid lines have slopes equal to —1.5. 
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Further, we fix mass and charge of dilaton black hole and change the mass of the SI scalar fields. We have noted the 
relation between mass m and the another pattern of decay; namely, the smaller m is the later begins to dominate this 
pattern. We remark that analogical relation was obtained from analytical considerations when the charge of dilaton 
black hole was assumed to be small. In our numerical simulations we get rid of this assumption. 



FIG. 4. Late-time behavior of the SI scalar field | ip \ max (maxima of the oscillations) at future timelike infinity for different 
masses m as a function of time t. The thin solid line has slope equals to —1.5, mass of dilaton black hole M = 0.5 and charge 



In our work we have investigated the intermediate asymptotic behavior of SI scalar field in the background of dilaton 
black hole being the spherical symmetric solution of the so-called low-energy string theory |lj] . As was revealed the 
late-time behavior of the charged massless perturbations depended on the spacetime parameters M, Q, while the 
intermediate asymptotic behavior of SI field depended only on the field's mass m. The numerical scheme has shown 
that at the very late-time the inverse power-law decay is replaced by another kind of the decay process, slower than 
any power-law. From our numerical calculations we have established that if one fixed the mass and charge of dilaton 
black hole M we get the relation: the smaller SI field mass is, the later begins the process. This very late-time 
behavior needs further analytic studies. The first step towards it was done in ]ll] , |T2| studying nearly extreme RN 
and Schwarzschild solutions. We have found that in the case of small Q comparing to M, our wave equation reduce 
to the form studied in Thus, it can be deduced that the asymptotic late-time behavior of SI scalar field is of 

the form of the oscillatory tail with the decay rate i -5 / 6 , which is identical with that in the nearly extreme RN or 
Schwarzschild background. 

Another problem for further investigations is the problem of late-time behavior of massless and SI scalar fields in 
the spacetime of extremal dilaton black hole. We hope to return to this problem elsewhere. 




Q = 0.45. 



IV. CONCLUSIONS 
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